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In order to test the role of large-scale magnetic fields in quasiperiodic oscillation phenomena ob¬ 
served in microquasars, we study oscillatory motion of charged particles in vicinity of a Schwarzschild 
black hole immersed into an external asymptotically uniform magnetic field. We determine the fun¬ 
damental frequencies of small harmonic oscillations of charged test particles around stable circular 
orbits in the equatorial plane of a magnetized black hole, and discuss the radial profiles of fre¬ 
quencies of the radial and latitudinal harmonic oscillations in dependence on the mass of the black 
hole and the strength of the magnetic field. We demonstrate that assuming relevance of resonant 
phenomena of the radial and latitudinal oscillations of charged particles at their frequency ratio 
3 : 2, the oscillatory frequencies of charged particles can be well related to the frequencies of the 
twin high-frequency quasi-periodic oscillations observed in the microquasars GRS 1915+105, XTE 
1550-564 and GRO 1655-40. 
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I. INTRODUCTION 

Investigation of the high-frequency quasi-periodic os¬ 
cillations (HF QPOs) observed in many black hole or 
neutron star low-mass X-ray binaries (LMXB) can open 
up prospects of understanding the phenomena occurring 
in the presence of strong gravity. Thanks to the pos¬ 
sibilities to achieve high precision in the measurements 
of the frequencies of observed oscillations one can get 
some useful information about the central object and the 
electromagnetic fields in its vicinity from the analysis of 
obtained frequencies. 

In the LMXB systems containing a black hole or a neu¬ 
tron star, the HF QPOs are sometimes observed in pairs 
of the upper and lower frequencies {v\j, of twin peaks 
in the Fourier power spectra. In so called microquasars, 
i.e., LMXB systems containing a black hole, the twin HF 
QPOs occurs at fixed frequencies that are usually at the 
exact 3 : 2 ratio [I]. The observed high frequencies are 
close to the orbital frequency of the marginally stable cir¬ 
cular orbit representing the inner edge of the Keplerian 
discs orbiting black holes, therefore, the strong gravity 
effects have to be relevant in explaining HF QPOs [2]. 
Usually, the Keplerian orbital and epicyclic (radial and 
latitudinal) frequencies of the circular geodesics of the 
Kerr geometry [3-6] are assumed in models explaining 
the HF QPOs in the black hole systems. Alternatively, 
the oscillations of tori [7] or tilted oscillating discs [8] are 
considered as explanation of the HF QPOs. The frequen¬ 
cies of the oscillating tori and tilted discs are related to 
the orbital and epicyclic geodesic frequencies, if the os¬ 
cillations are governed mainly by the gravity of the black 
hole [9]. 

It was suggested in [10] that resonances in oscillations 
of Keplerian accretion discs have to be relevant. This 
assumption seems to be confirmed due to the properties 


of the HF QPOs observed in the three Galactic micro¬ 
quasars GRS I9I5+105, XTE 1550-564, GRO 1655-40, 
where the 3 : 2 frequency ratio has been observed [2]. 
Note that there has been observed more complex fre¬ 
quency pattern in the microquasar GRS 1915+105, con¬ 
taining five frequencies, and more complex theory has to 
be involved in order to explain the whole frequency set 
[11, 12]. Here we focus our attention to the twin HF 
QPOs demonstrating the 3 : 2 frequency ratio. 

Unfortunately, the observed frequencies of the 3 : 2 HF 
QPOs in the three microquasars cannot be explained by 
a unique model based on the frequencies of the geodesic 
epicyclic motion, if the limits on the mass and spin of 
the black holes are taken into account [13, 14]. In the 
present paper, we devote our attention to the possibil¬ 
ity to explain the observed 3 : 2 HF QPOs in the three 
microquasars due to charged particles oscillating around 
stable circular geodesics in combined gravitational and 
electromagnetic fields of magnetized black holes. Obser¬ 
vational evidence for existence of a magnetic field around 
astrophysical black holes can be found in [15, 16]. 

The role of the magnetic field in the processes taking 
place in the black hole surroundings is relevant due to 
several reasons. Most of the observed black hole candi¬ 
dates have an accretion disc constituted from conduct¬ 
ing plasma which dynamics can generate a regular mag¬ 
netic field. Recently, it has been found that the center 
of the Galaxy has a strong magnetic field around a su- 
permassive black hole that is not related to an accretion 
disc [17]. Therefore, black holes can be also immersed 
in an external, large scale electromagnetic field that can 
have a complicated structure in vicinity of field source, 
but at large distance (asymptotically) in finite element of 
space, its character can be simple and close to a homo¬ 
geneous magnetic held. Moreover, it has been recently 
shown [18] that a black hole located near the equatorial 
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Figure 1: Quasi-periodic motion of a charged test particle around a Schwarzschild black hole immersed in an external uniform 
magnetic field (right) compared to the motion of an electrically neutral test particle moving along a circular orbit around the 
black hole (left). The charged particle starts its motion at the point on the circular geodesic of the uncharged particle. If 
the energy and angular momentum of the charged particle is slightly modified with respect to those of the circular geodesic, 
quasi-harmonic oscillations can occur. 


plane of a magnetar will be immersed in a nearly ho¬ 
mogeneous magnetic field if distance to the magnetar is 
large enough. Hereafter in this paper we will concentrate 
our attention on this particular and simplified case of a 
black hole immersed in an asymptotically uniform mag¬ 
netic field known as Wald solution for a magnetized black 
hole [19]. 

According to [20], where the magnetic coupling process 
has been considered, in the vicinity of a stellar mass, 
M ~ 10 Mq, or a supermassive, M ~ 10® Mq, black 
hole, the estimate of the strength of the magnetic field 
gives, B ^ 10® Gauss, or, B ~ 10"^ Gauss, respectively. It 
has to be stressed that in both cases the magnetic field 
does not change the gravitational background so that the 
metric tensor of the Schwarzschild black hole spacetime 
does not need a modification. In order to satisfy this test 
field approximation, the strength of the magnetic field in 
the vicinity of a black hole with the mass M has to fulfil 
the condition [21] 

S«-BG = g37jM0(^j~lO'»^Ga„ss. (1) 

The magnetic fields satisfying the condition of the test 
field approximation, B « Bq, have negligible effect on 
the motion of neutral particles. However, for the motion 
of charged test particles the influence of the magnetic 
field can be really large. For a charged test particle with 
charge q and mass m moving in vicinity of a black hole 
with mass M surrounded by an external asymptotically 
uniform magnetic field of the strength B, one can intro¬ 
duce a dimensionless quantity b that can be identified as 
relative Lorenz force [21] 

i, = !£!£^. (2) 


This quantity can be quite large even for weak magnetic 
fields due to the large value of the specific charge q/m [22, 
23] 


b 



B 


10® Gauss 


M \ 
10 Mq) ’ 


( 3 ) 

where nip is the proton mass and e is its charge. Accord¬ 
ing to the estimation (3) the influence of the magnetic 
field on the motion of charged particles cannot be ne¬ 
glected even for weak magnetic fields. 

Motion of charged test particles in combined gravita¬ 
tional and electromagnetic fields around black holes has 
been studied in a large variety of papers both for the 
Reissner-Nordstrom and Kerr-Newman solutions of the 
Einstein-Maxwell equations [5, 24-27] and black holes 
with external test electromagnetic fields [21, 23, 28-33]. 
Existence of the off-equatorial circular motion of charged 
particles in the combined electromagnetic and gravita¬ 
tional fields has been considered in [34-36]. Possible ex¬ 
istence of the off-equatorial levitating toroidal configura¬ 
tions has been also demonstrated in [37]. The astrophys- 
ically relevant energetic mechanisms such as the ’’mag¬ 
netic” Penrose process or the Blandford-Znajek mech¬ 
anism for black holes surrounded by a toroidal electric 
current have been studied in [38, 39]. 

The motion of charged test particles can be treated 
using Hamilton equations and related effective potential. 
For black holes with an external magnetic field, the mo¬ 
tion is generally of chaotic character, but in vicinity of 
the local minima of the effective potential the motion 
is simplified and corresponds to linear harmonic oscilla¬ 
tions or quasi-harmonic oscillation with frequencies close 
to those of the harmonic motion. The quasi-harmonic 
oscillations around a stable equilibrium location can be 
relevant for explanation of the HF QPOs observed in the 
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three microquasars GRS 1915+105, XTE 1550-564, and 
GRO 1655-40 [2, 40]. 

In this paper we consider motion of charged test parti¬ 
cles around a Schwarzschild black hole immersed in an ex¬ 
ternal uniform magnetic field. We look especially for ex¬ 
istence and properties of the harmonic or quasi-harmonic 
oscillations of charged particles in the magnetized black 
hole backgrounds, see Fig. 1. Many individual charged 
particles will probably not create complete accretion disc, 
since the fluid element of the accretion disk will be largely 
charge neutral, but the charged particle epicyclic motion 
can still be used for dynamics of ionized blob structures 
created by instabilities or by irradiation in otherwise neu¬ 
tral accretion disk [9]. 

Throughout the paper, we use the spacetime signature 
(—, +, +, +), and the system of geometric units in which 
G = 1 = c. However, for expressions having an astro- 
physical relevance we use the speed of light explicitly. 
Greek indices are taken to run from 0 to 3. 


The commuting Killing vector = djdij) generates ro¬ 
tations around the symmetry axis. Gonsequently, the 
only nonzero covariant component of the potential of the 
electromagnetic field takes the form [41] 

B B 

^ 500 = -2 sin^ e. (9) 

The charged test particle motion is described by the 
Lorentz equation 

TO—— = qFjJ'u^, (10) 

dr 

where is the four-velocity of the particle with the mass 
TO and charge q, normalized by the condition = — 1, 
T is the proper time of the particle, and 
is the antisymmetric tensor of the electromagnetic field. 

A. Hamiltonian and effective potential 


II. CHARGED PARTICLE DYNAMICS 


Using the general Hamiltonian formalism, we describe 
dynamics of a charged particle with charge q ^ 0 in 
the vicinity of the Schwarzschild black hole embedded 
in external asymptotically uniform magnetic field. The 
equations of motion of neutral particles can be obtained 
by taking q = 0. 

The line element of the Schwarzschild black hole space- 
time with mass M reads 


ds2 = -f{r)df + f-\r)dr^ + r^{d0^ + sin^ Odcj)'^), (4) 
where the function /(r) takes the form 


. , 2M 

/w = i-—■ 


( 5 ) 


Hereafter, we put M = 1, i.e., we use dimensionless radial 
coordinate r (and time coordinate t). 

The stationarity and axial symmetry of the 
Schwarzschild spacetime imply existence of the timelike 
and spacelike Killing vectors that satisfy the equation 


Ca;/3 + C/3;a = 0. (6) 

This implies that the solution for the four-vector poten¬ 
tial A^ representing the test electromagnetic field takes 
the following form [41] 

A'^ = Gie';) + G2ef,). (7) 

We consider the case of magnetic field which is uniform 
at the spatial infinity, having strength B there. The 
field is oriented perpendicularly to the equatorial plane of 
the black hole spacetime. Then the four-vector potential 
takes the form 

^" = 1 +)- ( 8 ) 


The Hamiltonian for the charged particle motion can 
be written in the form [41] 

Hp = - qAa){nfi - qA^) + i TO^ (II) 

where the kinematical four-momentum = mu^ is re¬ 
lated to the generalized (canonical) four-momentum 
by the relation 


that satisfy the Hamilton equations in the form 
dx'^ _ ^ _ dH d-K^ _ dH 
dC ^ ’ dC ' 


( 12 ) 


(13) 


The afhne parameter C, of the particle is related to its 
proper time r by the relation ^ = Tim. 

Due to the symmetries of the Schwarzschild spacetime 
(4) and the asymptotically uniform magnetic field (9), 
one can easily find the conserved quantities that are the 
energy and the axial angular momentum of the particle 
and can be expressed as 


E = -7rt = TO/(r)^, (14) 

L = = + (15) 


The dynamical equations for the charged particle mo¬ 
tion in the Cartesian coordinates can be found by the 
coordinate transformations 


X = r cos((/)) sin(d ), y = r sin((/)) sin(d ), z = r cos{9). 

( 16 ) 

Introducing for convenience the specific parameters, 
energy f, axial angular momentum £, and magnetic pa¬ 
rameter B, by the relations [21] 


£ = 


E 

to ’ 


C 


L 

to ’ 



(17) 
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Figure 2: Sections (y = 0 , z =const.) of the effective potential Vea{x,y,z;C,B) taken at the equatorial plane z = 0 (solid 
curves) and at the spatial infinity a —^ oo (dashed curves) for a charged test particle moving around a Schwarzschild black hole 
immersed in an uniform magnetic field. We compare uncharged case, B = 0 (left), with negative, B < 0 (middle), and positive, 
B > 0 (right), configurations. Values £ = 4 and = ± 0.005 are used for the graphs; the minima of VeH(x,z —>■ oo) (dotted 
vertical line) is at r = 20%/2. The influence of the external magnetic field becomes decisive for large values of coordinate x, 
for small values of x, the role of the magnetic field is suppressed by gravity in the equatorial plane and its vicinity, and by the 
centrifugal forces related to the angular momentum for z ^ oo, see eq. (19) 


one can rewrite the Hamiltonian (11) in the form 

1 11 Tn ^ 

+ ^pI + 2^^ [yesir,0) - £^] , (18) 

where Ves{r, 0; £, B) denotes the effective potential given 
by the relation 


charge of the particle is taken to be positive q > 0, the mi¬ 
nus configuration B < 0 corresponds to the vector of the 
magnetic field B pointing downwards, while plus config¬ 
uration B > 0 corresponds to the vector of the magnetic 
field B pointing upwards the z-axis. 

The charged particle motion is limited by the energetic 
boundaries given by 


Ves{r, 0) = /(r) 


1 + 


C 


rsin0 


— Br sin( 


(19) 


The terms in the parentheses corresponds to the cen¬ 
tral force potential given by the specific angular momen¬ 
tum C, and electromagnetic potential energy given by 
the magnetic parameter B. 

The effective potential (19) shows clear symmetry 
{C,B) O {—C,~B) that allows to distinguish the fol¬ 
lowing two situations 

- minus configuration, here £ > 0, H < 0 (equiva¬ 
lent to £ < 0,H > 0) - magnetic field and angular 
momentum parameters have opposite signs and the 
Lorentz force is attracting the charged particle to 
the z-axis, towards the black hole. 

-I- plus configuration, here £ > 0, H > 0 (equivalent 
to £ < 0,H < 0) - magnetic field and angular mo¬ 
mentum parameters have the same signs and the 
Lorentz force is repulsive, acting outward the black 
hole. 


The positive angular momentum of a particle £ > 0 
means that the particle is revolved in the counter¬ 
clockwise motion around the black hole, see Fig. 1. If 


£^ = V,s{r,e-,C,B). (20) 

Let us properly investigate the features of the effective 
potential (19) represented in Fig. 2 that enables us to 
demonstrate general properties of the charged particle 
dynamics, avoiding the necessity to solve the equations 
of motion [21]. 

The axial symmetry of the background of the combined 
gravitational and magnetic fields implies independence 
of the effective potential VeS on the coordinate cf which 
allows us to examine Ves{f,0) as a 2D function of vari¬ 
ables r, 9 or Cartesian x, z coordiantes (16). The effective 
potential is positive outside the black hole horizon, and 
diverges at the horizon r = 2. The region within the 
horizon and divergent points is excluded from our inves¬ 
tigation. 

The stationary points of the effective potential 
14ff(r, 0) function, where maxima or minima can exist, 
are given by the equations 

drVes{r,e-,C,B) = 0, deV,sir,0-,C,B)=O. (21) 

The second of the extrema equations (21) gives 9 = -k 12. 
In another words, all extrema of the 14ff(r, 0); £, B func¬ 
tion are located in the equatorial plane and there is no 
off-equatorial extremum for the charged particles orbiting 
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the Schwarzschild BH immersed in the uniform magnetic 
field, in contrast to the dipole magnetic field case [34]. 

The first extrema equation of (21) leads to a polyno¬ 
mial equation of the fifth order in the radial coordinate 

-3) + 2CBr'^ - - 1) - = 0, (22) 

which has generally five complex roots. Real roots of (22) 
above the black hole horizon determine maxima, minima 
and inflex points of the T4ff(r, 0 = 7r/2;£,S) function. 
Such extrema give stable (minima) and unstable (max¬ 
ima) equilibrium positions for the circular particle mo¬ 
tion, i.e. stable or unstable circular orbits. The inflex 
points give the marginally stable circular orbits. 

The equation (22) is quadratic with respect to the spe¬ 
cific angular momentum C and hence the circular orbits 
can be determined by the relation 

— Br^ '^tF 

C = C^±{r-B) = -—, (23) 

r — 6 

where we have introduced new function 


F{r-B) = A/BV2(r-2)2-Hr-3. (24) 

The positive branch of (23), the solution £e-i- > 0, can 
exist at the whole r > 2 region and determines both the 
stable and unstable circular orbits; for the plus configu¬ 
rations with B > Q, both functions /1 e± > 0 can exist in 
the 2 < r < 3 region, see Fig. 3. The negative branch of 
(23), the solution £e-j is responsible for maxima of the 
effective potential 14ff only. 

The local extrema of the /1 e±(^) function determine 
the innermost stable circular orbits (ISCO), their radius, 
angular momentum, and energy. They are determined 
by the condition 9^14®(r, 0; £, S) = 0 that implies the 
relation 


£E(ex)(r;S)=-2Sr + G, (25) 


where 


G{r;B) = \J B'^r'^ (5r^ — 4r -|- 4) -|- 2r. (26) 

The ISCO is located at the point of intersection of the 
functions /lE(ex) and £e+, see Fig. 3., which leads to the 
condition for the ISCO radius given in the form 

6-r+2B(r-6)G-2B^r (2r^ - 9r^ + 8r - 12) = 0. (27) 

Contrary to the uncharged test particle motion with 
B = 0, the motion of charged particles in the uniform 
magnetic field is always bounded in the r direction near 
the equatorial region 0 ~ 7r/2, because of the term B^r^ 
in effective potential (19), which unlimitedly grows with 
r —> oo. However, the energetic boundary (20) for 
the charged particle motion can be open in the polar 
direction of the black hole spacetime, the z direction, 
making the charged particles able to escape to infinity 


along the z axis. The energetic condition for the particles 
which can escapes to infinity reads 

c ^ c f 1 for H > 0 

t _ iflat(min) - | + 4g£ foj. g < Q 

where £flat(min) is the minimal energy of the charged par¬ 
ticle at infinity, see Fig. 4. 

If the energetic boundary (20) forms a closed curve, 
the motion of charged particles is trapped, and their tra¬ 
jectory is restricted to a toroidal-like region governed by 
the effective potential. The condition for formation of 
the closed “lake” like energy boundaries that are related 
to the local minima of the effective potential Ves{x,z) 
obviously reads 


^ < ^’flat(min)- (29) 

One can deduce that the trapped states of the oscillating 
charged particle can exist, if 

Cei<C< Ce2 (30) 


where the so called “lake” angular momentum functions 
£Li{r]B) and £l 2 (?';S), are the solutions of the equal¬ 
ity condition £ = i?flat(min)- For B > 0, the functions 
£Li{r;B) and £L 2 {r',B) take the form 


Br'^{r — 2) ± ry/2{r — 2) 

£li.L2 = -z- 

r — 2 

while for H < 0, they take the form 


—Br^{r -I- 2) ± r^j2{r — 2 4H^r^) 

£li,L2 = -z- 


(31) 


(32) 


The functions /iLi(r), £L 2 ('r) are represented in Fig. 3, 
for e = 0 and H = ± 0.1. 

The energy of the circular orbits of the charged parti¬ 
cles, given by the local minima and maxima of the effec¬ 
tive potential, is represented in Fig. 4. Typical sequences 
of the related charged particle trajectories, for given spe¬ 
cific angular momentum £, are represented in Figs. 5 
and 6. 


B. Curled trajectories 

Using the conservation of the angular momentum (15) 
of the motion of a charged particle around Schwarzschild 
black hole in uniform magnetic field, the equation of mo¬ 
tion (13) for the axial coordinate (j) can be written in the 
form 


^=^-B. (33) 

For uncharged particles, the case H = 0 (£ > 0, r > 0), 
the right side of the above equation is always positive and 
hence the coordinate tp can only increase. But for positive 
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Figure 3: Circular orbits of charged particles. The extrema function ZlE±(r) (thick) determines local extrema of the effective 
potential Feff(r, 0), i.e., the circular orbits in the equatorial plane. -Cisco governs the innermost stable circular orbit (relative to 
radial perturbations), the stable (unstable) circular orbits are located at r > nsco {r < risco)- The functions -CLi(r),-CL 2 (r) 
(dashed) govern the trapped states. The function €.E(ex){'f", B) (dotted) governs the extrema of the CE±{r) function. The 
negative branch of the extrema function, £.-E-{r) (thick dashed curve), and the curve (dot-dashed) governing the 

existence of the trajectories with curls, can both exist only for positive values of the magnetic parameter, B > 0. Regions 
where trapped states (bounded orbits only) can exist are shaded. All curves are given for neutral B — 0 (left), negative B < 0 
(middle) and positive B > 0 (right) values of the charged test particle magnetic parameter. Numbers denote the examples of 
different types of trajectories illustrated in Figs. 5 and 6. 




Figure 4: Energetics of the circular orbits. Local extrema of the effective potential V),ff (r, 8) determine the circular orbits in the 
equatorial plane and are given as functions of the specific angular momentum £. Thick solid curves correspond to the maxima 
(unstable orbits) and minima (stable orbits) of the effective potential Vef[(r,8). Minimum of the effective potential for the flat 
spacetime f’flat(min) (governing escape to infinity) is represented by the dashed line. For B > 0, the curly motion is allowed for 
particles with parameters above the function f’«(-C) (dot-dashed line). Regions where the trapped states can exist are shaded. 
Numbers denote the examples of different types of trajectories illustrated in Figs. 5 and 6. 


values of the electromagnetic interaction parameter, B > 
0, the coordinate (j) can decrease, if 

C> C4r;B) = Br^. (34) 

This condition can be expressed in terms of the energy 
condition [21] 


(35) 


Decrease of the coordinate (j) during the charged particle 
motion leads to the epicyclic motion with curled trajec¬ 
tories. The ’’curling” functions, £*(r) and £*{C), are 
illustrated in Figs. 3 and 4, respectively. The transition 
from a trajectory without curls to those demonstrating 
the curly character is illustrated for a constant specific 
angular momentum £ in Fig. 6. Note that the curled 
trajectories are the most illustrative example of the in- 


£ > £^{C;B) = \ll-2y/BjC. 
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Figure 5: List of charged test particle equatorial trajectories with specific angular momenta £ = 8 around a Schwarzschild 
black hole immersed in an uniform magnetic field - the magnetic parameter B = —0.1, see Fig. 3. (middle). The particle starts 
from different initial radii r clarifying different shapes of its trajectory. Dotted circle represents stable circular orbit for given 
specific angular momenta C. 



Figure 6 : Curling trajectories. List of charged test particle equatorial trajectories with specific angular momenta C = 4.5 
around a Schwarzschild black hole immersed in an uniform magnetic field - the parameter B — 0.1, see Fig. 3. (right). 


fluence of the magnetic field. 


C. Trajectories of charged particles 

Generally, the motion of a charged particle around a 
Schwarzschild black hole immersed in an external uni¬ 
form magnetic field (Wald solution) is chaotic [36, 42]. 
However, trajectories of charged particles, which are close 
to the minima of the effective potential corresponding to 
the stable circular orbits, still have a regular (harmonic) 
character [36]. Moreover, the trajectories entirely located 
in the equatorial plane are also regular, and the chaotic 


behaviour occurs due to varying the inclination angle 6*0 
from the equatorial plane [42]. Here we will focus on the 
regular (harmonic), or quasi-harmonic bounded orbits lo¬ 
cated close to the equatorial plane, with initial inclination 
angle 0o ~ 7 ’'/ 2 , see Fig. 1. 

Due to the behaviour of the effective potential 
Ves{x,z), one can distinguish four different types of the 
energetic boundary (20). The differing energetic bound¬ 
aries can be seen as dashed curves in the third line in 
Fig. 7. The first case corresponds to non existing inner 
and outer boundary - the particle can be captured by the 
black hole or escape to infinity, see Fig. 7(d). The second 
case corresponds to the situation with an outer boundary 
















(a) collapse 


(b) bounded 


(c) bounded with curls 


(d) escape 


Figure 7: Typical trajectories of charged test particles (solid curves) orbiting a Schwarzschild black hole (shaded circle) immersed 
in an uniform magnetic field. The value of the electromagnetic interaction parameter is chosen as B = ±0.1. We give x-y and 
x-z sections (first and second line) of complete 3D x-y-z particle trajectory (fourth line). Due to the conservation of the particle 
specific energy £ and specific angular momentum C, the trajectory in 4D configuration space {t, x, y, z) can be represented in 
2D x-z graph (third line), where we also plotted the boundary of the particle motion given by the effective potential (dashed 
curve). 


- the charged particle must be captured by the black hole, 
see Fig. 7(a). The third case corresponds to the situation 
when both inner and outer boundaries exist - the charged 
particle is trapped in some region forming a toroidal re¬ 
gion around the black hole horizon, corresponding to a 
potential “lake” around the black hole, see Figs. 7(b,c). 


The fourth case corresponds to the existence of the inner 
boundary - the particle cannot fall into the black hole, 
but it must escape to infinity. 

We can also distinguish four different types of the 
charged particle regular motion. We start with the col¬ 
lapse trajectories which end at the horizon (we discuss 
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motion only above horizon, at r > 2). Undoubtedly, the 
collapsed trajectories can appear only in the first and sec¬ 
ond cases of the energetic boundary that are to enable 
particles to reach the horizon. An example of such col¬ 
lapsed trajectory is represented in Fig. 7(a). Then there 
are escaped trajectories, which end at the spatial infin¬ 
ity, z ^ oo, evolving along the magnetic field lines. Es¬ 
caped trajectories can appear only in the first and fourth 
cases of the energetic boundary, enabling the possibility 
to escape the strong gravity black hole region along the 
z-axis. An example of such trajectory can be found in 
Fig. 7(d). The escape along the z-axis appears to be 
uniquely related to the charged particle motion in the 
uniform magnetic field. 

Recall that motion of electrically neutral test particle 
around a Schwarzschild black hole is always restricted to 
a central plane, and freely moving particle can escape 
to infinity only in this plane a; —>■ oo. The last type of 
trajectories are the bounded trajectories -- their radial 
coordinate is bounded by apoapsis and periapsis rp, 
such as Tp < r < Ta. There can be two types of 
the bounded trajectories, namely those with and with¬ 
out curls - examples of such trajectories can be found 
in Figs. 7(c,b). The trajectories with curls are specific 
to the charged particle motion in the asymptotically uni¬ 
form magnetic fields, and can not occur in the motion of 
uncharged particles. A special type of the bounded or¬ 
bits are the stable circular orbits. We demonstrate that 
the ISCO of the charged particle circular orbits is always 
located inside the ISCO of neutral particles, located at 
r = 6. 


III. HARMONIC OSCILLATIONS 



Figure 8: Locally measured radial (horizontal) a;r, latitudinal 
(vertical) uje, Keplerian and Larmor wl angular frequen¬ 
cies for charged particle moving in vicinity of a stable circular 
orbit in the gravitational Schwarzschild field combined with 
the uniform magnetic field. Note that the Larmor frequency 
is related exclusively to the magnetic field and is thus relevant 
at large distance from the black hole. 


and electromagnetic background are then given by 


-2r {Br^ - Cf - 2r^] (39) 

where C = £e-i- is the specific angular momentum at the 
circular orbit (23). 

Recall that there exist the third fundamental angular 
frequency of the epicyclic particle motion, namely the 
Keplerian (axial) frequency given by 


If a charged test particle is slightly displaced from the 
equilibrium position located in a minimum of the effec¬ 
tive potential I4ff (u ^) at tq and = 7r/2, corresponding 
to a stable circular orbit, the particle will start to oscil¬ 
late around the minimum realizing thus epicyclic motion 
governed by linear harmonic oscillations. For harmonic 
oscillations around the minima of the effective poten¬ 
tial I4ff, the evolution of the displacement coordinates 
r = tq + 5r,9 = 9 q + 59 is governed by the equations 

5r -\- 5r = 0, 59 + ojg 59 = 0, (36) 


where dot b = do/dr denotes derivative with respect to 
the proper time r of the particle, and locally measured 
angular frequencies of the harmonic oscillatory motion 
are given by [41] 






eff 


dr"^ 


ujj = 


1 


eff 


r‘^f{r) 89^ 


(37) 


For an alternative definition of the epicyclic harmonic 
motion see [29]. 

Locally measured latitudinal (vertical) ojg and radial 
(horizontal) oj,- angular frequency of the harmonic oscil¬ 
lations of charged particles in the combined gravitational 


uj^ = ^ = U^ = —-B. 
dr 


(40) 


With the uniform magnetic field itself the so called Lar¬ 
mor angular frequency wp is associated, being given by 
the relation 


WL = — = 2B. (41) 

m 

Obviously, the Larmor angular frequency wp is not de¬ 
pendent on the radial coordinate r and it is fully relevant 
in large distances from the black hole where the uniform 
magnetic field becomes to be crucial. For illustration of 
definition of the frequencies see Fig. 8. 

Behaviour of the locally measured angular frequencies 
ujr{r), U!g{r), u!^{r), and wp(r), as functions of the radial 
coordinate r, is demonstrated for significant values of 
the electromagnetic interaction parameter, B = 0 and 
B = ± 0.01,0.1, in Fig. 9. For small radii, r > nsco, 
we see strong gravitational influence on the angular fre¬ 
quencies (except the Larmor frequency, of course), for 
large radii r ^ Usco fhe influence of the uniform mag¬ 
netic field is prevailing. Close to the ISCO, the latitu¬ 
dinal uig frequency is always larger than the radial Wr 
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r 


r 


Figure 9: Radial profiles of the locally measured fundamental frequencies cur, u)e, and cjl for the charged particle motion. 
We compare the case of uncharged particle motion with S — 0 (left) to the charged particle motion with B < 0 (middle) and 
B > 0 (right). We give the radial profiles for B — i: 0.01 and B = ± 0.1; the same behaviour occurs for other values of the 
parameter B. 


frequency, ujg > uJi, but since Wr wl and we —>■ 0 
as coordinate r is increasing, there always exists radius 
ri:i where Wr = uig. The radial prohles of the latitudinal 
iog and radial Wr angular frequencies do not demonstrate 
significant dependence on the sign of the parameter B. 
The radial angular frequency Wr is approaching the Lar- 
mor frequency wl, from above for B < 0 and from below 
for S > 0, as the coordinate r is increasing Wr —?► wl. 
The Keplerian frequency is approaching asymptoti¬ 
cally the Larmor frequency wl from above for B < 0, but 
for > 0 it vanishes asymptotically. 

Behaviour of the fundamental frequencies uj,-,u}g,uj^ 
and their ratios can help us to distinguish different shapes 
of charged particle epicyclic orbits in the vicinity of a sta¬ 
ble circular orbit. According to the Newtonian theory of 
gravitation, all the frequencies are equal, = ojg = 
giving ellipse as the only possible bounded trajectory 
of a test particle around a gravitating spherically sym¬ 
metric body. For uncharged particles moving around a 


Schwarzschild black hole the relation < ujg = holds, 
and there exist a periapsis shift for bounded elliptic-like 
trajectory implying the effect of relativistic precession 
that increases with decreasing radius of the orbit as the 
strong gravity region is entered [3]. For charged parti¬ 
cles orbiting a Schwarzschild black hole immersed in an 
uniform magnetic field, the relativistic precession effect 
is modified by the presence of the magnetic field, and we 
can observe for the positive electromagnetic interaction 
parameters, S > 0, a special family of bounded orbits 
with curls. 

Spiral orbits which are resembling a toroidal solenoid, 
have to satisfy the condition (wr ~ wg) ^ w^, but 

this condition is not valid for charged particles orbiting a 
charged source in the weak gravity limit [43], and we have 
shown that it is not possible even for the charged particles 
orbiting a nonrotating Schwarzschild black hole placed 
in an uniform magnetic held. On the other hand, in the 
held of rotating Kerr black holes and naked singularities. 
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Figure 10: Radial profiles of the fundamental radial and latitudinal harmonic frequencies Vr (thin) and ug (thick) related to 
the static distant observers. Uncharged particle frequencies (left), are compared to the charged particle frequencies for B < Q 
(middle), and B > 0 (right). The gravitational attractor is a Schwarzschild black hole with M = 10 Mq. The latitudinal 
frequency radial profiles, ve{r), are restricted to the region of existence of the stable circular orbits, limited thus by the ISCO 
where the radial frequency v^(r) vanishes. 


the spiral orbits can exist because of the existence of 
relativistic orbits with low (Keplerian) angular velocity 
relative to distant static observers [27, 44]. 


A. Frequencies measured by distant observers 

The locally measured angular frequencies and 

given by 


around a Schwarzschild black hole immersed in an uni¬ 
form magnetic field given by 

02/ m r-6 + 4eVi7 + 8S3(r-2)r3F , 

= - 4 /. , .02 2^ -^: (45) 


r4 (1 -hdBV) 




(46) 


_ 

’ r3(2BV2(r-2)-h2BrF-hl)’ 


(47) 



(42) 


where we have used the function F{r\B) (24) and intro¬ 
duced the function 


where P G {r, 0, (/>}, are connected to the angular frequen¬ 
cies measured by the static distant observers, 17, by the 
gravitational redshift transformation 


0^ 


dop dr ojp 

dt dt f{r)£{r)' 


(43) 


where {f{r)£{r))~^ is the redshift coefficient, given by 
the function /(r) and the particle specific energy at the 
circular orbit £(r). 

If the fundamental frequencies of the small harmonic 
oscillations related to the distant observers, £lp, are ex¬ 
pressed in the physical units, their dimensionless form 
has to be extended by the factor c^/GM. Then the fre¬ 
quencies of the charged particle radial and latitudinal 
harmonic oscillations measured by the distant observers 
are given by 


^/3 


1 


27r GM 


■ £if. 


(44) 


The dimensionless radial Hr, latitudinal and Kep¬ 
lerian angular frequencies related to the distant ob¬ 
servers are for the charged particle harmonic oscillations 


H{r-B) = 2SV3 - -t SeV + r-A. (48) 

In the models of the HF QPOs in the black hole systems, 
the epicyclic oscillatory frequencies are usually rel¬ 

evant [2]. Therefore, we focus our attention here on the 
frequencies of the radial and latitudinal harmonic oscil¬ 
lations. Frequencies i/i.(r),//^(r) related to the distant 
observers, as functions of the radial coordinate r, are 
plotted in Figs. 10 and 11. 

In the pure Schwarzschild black hole spacetimes {B = 
0), the harmonic oscillations have frequencies related to 
static distant observers given by expressions that are rel¬ 
atively very simple 

^?(geo)(^) = ^liseo)ir) = (49) 

It is quite interesting that the latitudinal frequency of 
the epicyclic geodetical motion equals to the 

charged particle latitudinal frequency in the combined 
Schwarzschild and uniform magnetic fields, see eq. (46). 
This effect is caused by the direction of the Lorentz force 
acting on the charged particle (resulting from the ex¬ 
ternal uniform magnetic field) - in the equatorial plane 
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Figure 11: Frequencies r'6i,t'r of small harmonic oscillations measured by static distant observers are given for the charged 
particles in the combined Schwarzschild and uniform magnetic fields. Position of the resonant radii are also presented. 


where the oscillations occur, it acts only in the radial 
direction. The motion is free, influenced only by grav¬ 
itational and inertial forces in the vertical (latitudinal) 
direction. 


B. Properties of the radial profiles of frequencies of 
the charged particle harmonic oscillations 


Both cases show similar behaviour for small values of the 
B parameter, for larger values of the B parameter the 
important radii tend to the value rg —> 2 for S —^ oo, 
while for B —>■ —oo we obtain the critical radius 

ro ^ i (5 -b = 4.30278, (53) 

see Fig. 12. 


We discuss properties of the radial profiles of both ra¬ 
dial flj. and latitudinal (vertical) fig angular frequencies 
of the harmonic oscillations, they are relevant in the mod¬ 
els of the HF QPOs [2]. The behaviour of latitudinal har¬ 
monic angular frequency, Vt^{r\B) (46), is quite simple, 
being of the standard Newtonian character, and does not 
depend on the parameter B. 

The zero points of the radial profile of angular fre¬ 
quency of the radial harmonic oscillations are given by 
the condition 


O2(r;,B) = 0 (50) 

that corresponds to the limit on the existence of the sta¬ 
ble orbits, the ISCO radii (27). Below ISCO, the particles 
start to collapse into black hole horizon, and no radial 
oscillations can occur. The local extrema of the radial 
profiles of angular frequency of the radial oscillations are 
given by the condition 

(51) 

dr 

The coincidence condition for the radial and latitudinal 
angular frequencies 

n,{r-B)=ne{r-B) (52) 

implies the existence of only one coincidence radius ri,! 
for all values of the parameter B. 

Summary of the properties of the function can 

be found in Fig. 12 for both B < 0 and B > 0 cases. 


Source GRO 1655-40 XTE 1550-564 GRS 1915-bl05 


lAj [Hz] 447 — 453 273 — 279 

iZL [Hz] 295 — 305 179 — 189 

M/Mq 6.03 — 6.57 8.5 — 9.7 

a 0.65 — 0.75 0.29 — 0.52 


165 — 171 
108 — 118 
9.6 — 18.4 
0.98 — 1 


Table I: Observed twin HF QPO data for the three micro¬ 
quasars, and the restrictions on mass and spin of the black 
holes located in them, based on measurements independent of 
the HF QPO measurements given by the spectral continuum 
fitting [45, 46]. 


IV. CHARGED PARTICLE OSCILLATIONS AT 
RESONANCE 3:2 AND 2:3 RADII MODELLING 
THE TWIN HF QPOS IN THE MICROQUASARS 

The quasi-harmonic character of the motion of charged 
particles trapped in a toroidal space around the equato¬ 
rial plane of magnetized black holes suggests an inter¬ 
esting astrophysical application related to the HF QPOs 
observed in the LMXB systems containing a black hole 
or a neutron star, or in an active galactic nuclei. Some 
of the HF QPOs come in pairs of the upper and lower 
frequencies (z/y, z^l) of twin peaks in the Fourier power 
spectra. Since the peaks of high frequencies are close to 
the orbital frequency of the marginally stable circular or¬ 
bit representing the inner edge of Keplerian discs orbiting 
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Figure 12: Properties of the radial profile of the fundamental angular frequency of the radial harmonic oscillations, flr(r;Z3), 
of charged particles oscillating in the combined Schwarzschild and uniform magnetic fields, (left) Position of the ISCO in 
dependence on the parameter B that reveals from which radius the bounded motion of charged particles can start. The radial 
harmonic frequency vanishes at the ISCO. (middle and right) Regions where the particle motion is bounded and the expressions 
for the harmonic oscillations are relevant are greyed, region below the black hole horizon is hatched. Solid thick curve represents 
the maximum of the flr(r) function, various resonant radii (where fig •. = n : m) are also plotted in dependence on the 

parameter B. 


black holes (or neutron stars), the strong gravity effects 
must be relevant in explaining HF QPOs [2]. 


Before the twin peak HF QPOs have been discovered 
in microquasars, first by [47], and the 3 : 2 ratio has been 
pointed out, [10] suggested that these QPOs should have 
rational ratios, because of the resonances in oscillations 
of nearly Keplerian accretion disks; see also [29]. It seems 
that the resonance hypothesis is now well supported by 
observations, and that the 3:2 ratio {2v\j = 3r'L) is seen 
most often in the twin HF QPOs in the LMXB contain¬ 
ing black holes (microquasars). Unfortunately, neither of 
the recently discussed models based on the frequencies of 
the harmonic geodesic epicyclic motion (based on Kep¬ 
lerian discs or slender tori [2, 48]) is able to explain the 
HF QPOs in all the three microquasars assuming their 
central attractor to be a black hole [13]. 

We can assume applicability of the parametric reso¬ 
nance, discussed in [49], focusing attention to the case of 
the frequency ratios izg : = 3 : 2 or zzg : t'r = 2 : 3, as 

the observed values of the twin HF QPO frequencies for 
GRO 1655-40, XTE 1550-564 and GRS 1915-bl05 sources 
show clear ratio 

= 3 : 2 (54) 

for the upper v\j and lower frequencies, see Tab. I. 
We identify directly the frequencies v\j , z/p with vg , or 
zzr,z/g frequencies. In contrast to the standard resonance 
epicyclic model (without external magnetic field), the os¬ 
cillating charged particles in the combined Schwarzschild 


and uniform magnetic fields allow both frequency ratios 

zzg : zZf = 3 : 2, vg : = 2 : 3. (55) 

Since < r 2-.3 (see Figs. 11 and 12), we call the first 
resonance radius, where vg : v-^ = 3 : 2^ the inner one, and 
the second resonance radius, where ug : = 2 : 3, the 

outer one. Note that for the oscillating charged particles 
also the 1:1,1:2,2:1, or other, resonant frequency ra¬ 
tios can enter the play - these can be relevant in other 
sources where such frequency ratios are observed, see, 
e.g., [50], or for the other twin frequencies observed in 
the microquasar GRS 1915-1-105. 

The procedure of fitting the charged particle oscilla¬ 
tion frequencies to the observed frequencies is presented 
in Fig. 13, for all the three microquasars GRS 1915-1-105, 
XTE 1550-564, and GRO 1655-40. From the restric¬ 
tions on the spacetime mass parameter M for each of 
the sources, see Tab. L, we obtain simultaneously re¬ 
strictions on the external uniform magnetic field through 
restrictions on the magnetic parameter B. So far we have 
not considered the sources rotation that can modify the 
radial profiles of the vertical and horizontal frequencies - 
the effect of the black hole rotation is going to be stud¬ 
ied in a future work, but our preliminary results show 
that the fitting can be done even for the extremely fast 
rotating microquasar GRS 1915-1-105. 

We use the fitting procedure in the case of relatively 
slowly (a ~ 0.4) rotating XTE 1550-564 source as an 
useful example. Detailed results for the XTE 1550-564 
source are given in Tab. H and HI, where we present typ¬ 
ical values of the magnetic field related to astrophysically 
realistic situations. Strong magnetic fields are not neces¬ 
sary - for electrons, the required magnetic field strength 
Be- ^ 0.1 mGs is comparable to the magnetic field 
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B 

Be- [mGs] Bp+ [Gs] Rfb [Gs] 

- 3:2 0.063 — 0.110 

- 2:3 0.113 — 0.370 

-t 3:2 0.041 — 0.056 

-1 2:3 0.048 — 0.058 

0.17 — 0.26 0.31 — 0.47 17.38 — 26.60 

0.31 — 0.88 0.56 — 1.61 31.18 — 89.47 

0.11 — 0.13 0.20 — 0.24 11.31 — 13.54 

0.13 — 0.14 0.24 — 0.25 13.25 — 14.03 


Table II: Boundary values of the magnetic field parameter B estimated for the microquasar XTE 1550-564. First part of this 
table contains restrictions on the B parameter for slowly rotating source XTE 1550-564 given by the resonance oscillation model 
of twin HE OPOs, see Fig 13. The mass of XTE 1550-564 source implies limits on the magnetic field strength. We present such 
limits for electron Be-, proton Bp+ and partially ionized (one electron lost) iron atom Byb- in Gauss units in second part. 


B 

Q X 10® [esu] m X 10 ^® [g] n x 10® 

- 3:2 0.063 — 0.110 

- 2:3 0.113 — 0.370 

+ 3:2 0.041 — 0.056 

+ 2:3 0.048 — 0.058 

0.90 — 1.38 5.33 — 3.49 3.18 — 2.08 

1.61 — 4.63 2.97 — 1.04 1.78 — 6.19 

0.59 — 0.70 8.19 — 6.85 4.89 — 4.09 

0.69 — 0.73 7.00 — 6.61 4.18 — 3.95 


Table III: Boundary values of the magnetic field parameter B estimated for the microquasar XTE 1550-564. Last three columns 
represent the inverse estimations of the specific charge Q = q/m, mass m and number of particles n of the oscillating test object 
(or blob), for the fixed value of the strength of the magnetic field B = 10® Gs. We assume the one-electron ionization, in the 
other words the blob has to be slightly ionized. More stronger ionization increases the mass m and the number of particles in 
the blob for the corresponding value defining by the number of electrons removed away from initially neutral blob. The blob 
mass is expressed in units of the mass of the neutron n in the last column. 


strength in the heliosphere, for protons, Bp+ ~ 0.2 Gs 
it is comparable to Earth’s magnetic field at its surface, 
and for partially ionized (one electron lost) iron atom 
B-pe ~ 10 Gs is comparable to the magnetic field strength 
in Earth’s core [51]. 

V. CONCLUSIONS 

External uniform magnetic field B can strongly influ¬ 
ence the charged particle motion around a magnetized 
black hole. In dependence on initial inclination angle 
9q, the charged particle can escape to infinity along the 
z axis related to the direction of the magnetic field, it 
can start to move chaotically, or it can start to oscil¬ 
late in a harmonic or quasi-harmonic regime. Even small 
magnetic field, ^ 10 Gs, can significantly influence the 
charged particle trajectory, if the specific charge of test 
particle is large enough. The oscillatory motion can be 
curled, giving a clear resemblance to the Larmor preces¬ 
sion in a pure magnetic field. 

We have shown that assuming relevance of the reso¬ 
nant phenomena of the radial and vertical (latitudinal) 
oscillations at their frequency ratio 3:2, the oscillatory 
frequencies of charged particle in uniform magnetic field 
of a magnetized black hole can be well related to the fre¬ 
quencies of the twin 3 : 2 HE QPOs observed in the mi¬ 
croquasars GRS 1915-bl05, XTE 1550-564, GRO 1655- 
40. We can conclude that simple radial and latitudinal 
charged particle quasi-harmonic oscillations can be con¬ 
sidered as one of the possible explanations of the HE 
QPOs occurring in the field of compact objects. 

The results of the procedure of fitting the frequencies 


of the twin HE QPOs in the three microquasars for the 
charged particle oscillations are similar to those obtained 
in [40] in the framework of the string loop oscillation 
model [52], indicating some similarities in the dynamics 
of charged particles and string loops. 

The effect of even weak external uniform magnetic 
field of a magnetized Schwarzschild black hole on the 
radial and latitudinal frequencies of harmonic oscilla¬ 
tions around a stable equatorial orbit is substantial; the 
magnetic field influence can erase all “strange” effects, 
like those coming from the external dimension in the 
braneworld model [4, 5]. 

The single isolated charged particle dynamics can be 
imprinted in description of diluted collisionless plasma 
toroidal structures [37], the kinetic description of colli¬ 
sionless plasma can effectively reflect the epicyclic mo¬ 
tion of charged particles [53]. Full magnetohydrodynam¬ 
ics simulations at the present state are not able to reflect 
the high-frequency QPOs, only low-frequency QPOs have 
been reflected in some models [54]. The charged particle 
epicyclic motion can still be used for dynamics of ionized 
blob structures created by instabilities or by irradiation 
in otherwise neutral accretion disk [9]. 

We can conclude that the charged particle quasi¬ 
harmonic oscillation model has to be considered seriously, 
but deeper details have to be studied and the results have 
to be confronted with the results of modelling other sig¬ 
natures of the strong gravity related to the optical phe¬ 
nomena, such as the spectral line profiles or light curves. 
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Figure 13: The upper oscillation frequency v\j at the 3:2 or 2:3 resonance radii, calculated for both the negative B < 0 (—) 
and positive 13 > 0 (+) values of the magnetic parameter, and compared to the mass-limits obtained from observations of the 
slowly rotating microquasars XTE 1550-564. The other two, fast rotating, microquasars GRO 1655-40 and GRG 1915-105 are 
also presented. 


Acknowledgments 

The authors would like to express their acknowledg¬ 
ments for the Institutional support of the Centre for The¬ 
oretical Physics and Astrophysics at the Faculty of Phi¬ 


losophy and Science of the Silesian University in Opava. 
The authors acknowledge the Albert Einstein Centre for 
Gravitation and Astrophysics supported by the Czech 
Science Foundation grant No. 14-37086G and the internal 
student grant of the Silesian University SGS/23/2013. 


[1] J. E. McGlintock, R. Narayan, S. W. Davis, L. Gou, 
A. Kulkarni, J. A. Orosz, R. F. Penna, R. A. Remil- 
lard, and J. F. Steiner, Glassical and Quantum Gravity 
28, 114009 (2011), 1101.0811. 

[2] G. Torok, M. A. Abramowicz, W. Kluzniak, and 
Z. Stuchh'k, Astron. Astrophys. 436, 1 (2005). 

[3] L. Stella and M. Vietri, Physical Review Letters 82, 17 
(1999), astro-ph/9812124. 

[4] A. Kotrlova, Z. Stuchh'k, and G. Torok, Glassical and 
Quantum Gravity 25, 225016 (2008), 0812.0720. 

[5] Z. Stuchhk and A. Kotrlova, General Relativity and 
Gravitation 41, 1305 (2009), 0812.5066. 

[6] J. Schee and Z. Stuchh'k, General Relativity and Gravi¬ 
tation 41, 1795 (2009), 0812.3017. 

[7] L. Rezzolla, S. Yoshida, and O. Zanotti, Mon. Not. R. 
Astron Soc. 344, 978 (2003), astro-ph/0307488. 

[8] S. Kato, Publications of the Astronomical Society of 
Japan 60, 111 (2008), 0709.2467. 

[9] Z. Stuchh'k, A. Kotrlova, and G. Torok, Astron. Astro¬ 
phys. 552, AlO (2013), 1305.3552. 

[10] W. Kluzniak and M. A. Abramowicz, ArXiv Astrophysics 
e-prints (2001), astro-ph/0105057. 

[11] Z. Stuchhk, P. Slany, G. Torok, and M. A. Abramowicz, 
Phys. Rev. D 71, 024037 (2005), gr-qc/0411091. 

[12] Z. Stuchh'k, P. Slany, and G. Torok, Astron. Astrophys. 
463, 807 (2007). 

[13] G. Torok, A. Kotrlova, E. Sramkova, and Z. Stuchh'k, 
Astron. Astrophys. 531, A59 (2011), 1103.2438. 

[14] G. Bambi, Journal of Gosmology and Astroparticle 
Physics 9, 014 (2012), 1205.6348. 


[15] M. de Kool, G. V. Bicknell, and Z. Kuncic, Publications 
Astronomical Society of Australia 16, 225 (1999). 

[16] J. M. Miller, J. Raymond, A. Fabian, D. Steeghs, 
J. Homan, G. Reynolds, M. van der Klis, and R. Wij- 
nands. Nature 441, 953 (2006), astro-ph/0605390. 

[17] R. P. Eatough and et ah. Nature 501, 391 (2013), 
1308.3147. 

[18] J. Kovaf, P. Slany, G. Gremaschini, Z. Stuchhk, 
V. Karas, and A. Trova, Phys. Rev. D 90, 044029 (2014), 
1409.0418. 

[19] R. M. Wald, Phys. Rev. D 10, 1680 (1974). 

[20] M. Y. Piotrovich, N. A. Silant’ev, Y. N. Gnedin, and 
T. M. Natsvlishvili, ArXiv e-prints (2010), 1002.4948. 

[21] V. P. Frolov and A. A. Shoom, Phys. Rev. D 82, 084034 
(2010), 1008.2985. 

[22] M. Y. Piotrovich, N. A. Silant’ev, Y. N. Gnedin, and 
T. M. Natsvlishvili, Astrophysical Bulletin 66, 320 
( 2011 ). 

[23] V. P. Frolov, Phys. Rev. D 85, 024020 (2012), 1110.6274. 

[24] R. Rufhni, in Black Holes (Les Astres Occlus), edited by 
C. Dewitt and B. S. Dewitt (1973), pp. 451-546. 

[25] D. Pugliese, H. Quevedo, and R. Rufhni, Phys. Rev. D 
83, 104052 (2011), 1103.1807. 

[26] J. Bicak, Z. Stuchhk, and V. Balek, Bulletin of the As¬ 
tronomical Institutes of Gzechoslovakia 40, 65 (1989). 

[27] V. Balek, J. Bicak, and Z. Stuchhk, Bulletin of the As¬ 
tronomical Institutes of Gzechoslovakia 40, 133 (1989). 

[28] A. R. Prasanna, Nuovo Cimento Rivista Serie 3, 1 (1980). 

[29] A. N. Aliev and D. V. Galtsov, General Relativity and 
Gravitation 13, 899 (1981). 
































16 


[30] A. M. A1 Zahrani, V. P. Frolov, and A. A. Shoom, Phys. 
Rev. D 87, 084043 (2013), 1301.4633. 

[31] G. Preti, Phys. Rev. D 70, 024012 (2004). 

[32] A. A. Abdujabbarov, B. J. Ahmedov, and N. B. Jurayeva, 
Phys. Rev. D 87, 064042 (2013). 

[33] P. Bakala, E. Sramkova, Z. Stuchh'k, and G. Torok, Glas- 
sical and Quantum Gravity 27, 045001 (2010). 

[34] J. Kovaf, Z. Stuchhk, and V. Karas, Classical and Quan¬ 
tum Gravity 25, 095011 (2008), 0803.3155. 

[35] J. Kovaf, O. Kopacek, V. Karas, and Z. Stuchh'k, Classi¬ 
cal and Quantum Gravity 27, 135006 (2010), 1005.3270. 

[36] O. Kopacek, V. Karas, J. Kovaf, and Z. Stuchh'k, The 
Astrophysical Journal 722, 1240 (2010), 1008.4650. 

[37] C. Cremaschini, J. Kovaf, P. Slany, Z. Stuchhk, and 
V. Karas, The Astrophysical Journal Supplement 209, 
15 (2013), 1309.3979. 

[38] S. M. Wagh and N. Dadhich, Phys. Rep. 183, 137 (1989). 

[39] L.-X. Li, Phys. Rev. D 61, 084016 (2000), astro- 
ph/9902352. 

[40] Z. Stuchhk and M. Kolos, Phys. Rev. D 89, 065007 
(2014), 1403.2748. 

[41] R. M. Wald, General relativity (University of Chicago 
Press, Chicago, 1984). 

[42] O. Kopacek and V. Karas, The Astrophysical Journal 
787, 117 (2014), 1404.5495. 

[43] J. Kovaf, European Physical Journal Plus 128, 142 
(2013). 


[44] Z. Stuchh'k, Bulletin of the Astronomical Institutes of 
Czechoslovakia 31, 129 (1980). 

[45] R. Shafee, J. E. McClintock, R. Narayan, S. W. Davis, L.- 
X. Li, and R. A. Remillard, Astrophysical Journal Letters 
636, L113 (2006), astro-ph/0508302. 

[46] R. A. Remillard and J. E. McClintock, Annual Review 
of Astronomy and Astrophysics 44, 49 (2006), astro- 
ph/0606352. 

[47] T. E. Strohmayer, ArXiv Astrophysics e-prints (2001), 
astro-ph/0104487. 

[48] L. Rezzolla, S. Yoshida, T. J. Maccarone, and O. Zan- 
otti, Mon. Not. R. Astron Soc. 344, L37 (2003), astro- 
ph/0307487. 

[49] L. D. Landau and E. M. Lifshitz, Mechanics (Oxford: 
Pergamon Press, 1969). 

[50] P. Lachowicz, B. Czerny, and M. A. Abramowicz, ArXiv 
Astrophysics e-prints (2006), astro-ph/0607594. 

[51] Y. B. Zeldovich, A. A. Ruzmaikin, and D. D. Sokoloff, 
Maqnetic Fields in Astrophysics (Gordon & Breach, New 
York, 1983). 

[52] M. Kolos and Z. Stuchhk, Phys. Rev. D 88, 065004 
(2013), 1309.7357. 

[53] C. Cremaschini and Z. Stuchhk, Phys. Rev. E 87, 043113 
(2013). 

[54] M. Machida and R. Matsumoto, Publications of the As¬ 
tronomical Society of Japan 60, 613 (2008), 0802.3779. 


